The large-scale structure of passive scalar turbulence 
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We investigate the large-scale statistics of a passive scalar transported by a turbulent velocity 
field. At scales larger than the characteristic lengthscale of scalar injection, yet smaller than the 
correlation length of the velocity, the advected field displays persistent long-range correlations due 
to the underlying turbulent velocity. These induce significant deviations from equilibrium statistics 
for high-order scalar correlations, despite the absence of scalar flux. 
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Turbulent flows are systems far away from equilib- 
rium, characterized by a flux of energy through a range 
of scales. This cascade process is often accompanied 
by strongly non-gaussian statistics and nontrivial scaling 
properties 0. Passive scalar turbulence is no exception, 
in this respect. Its understanding has recently undergone 
remarkable progress, with results thoroughly reviewed in 
Ref. |2(. A passive scalar field 9(x,t), like dilute dye 
concentration or temperature in appropriate conditions, 
transported by an incompressible velocity field v(x,t), 
evolves according to 



v-ve = kA6 + f , 



(i) 



where / is a source of scalar fluctuations that acts at a 
lengthscale L. The velocity field is assumed to be tur- 
bulent and characterized by a self-similar statistics (e.g. 
(v(x + r) — v{x)) ~ r 1 / 3 according to Kolmogorov's 1941 
theory) in the range of scales delimited above by the ve- 
locity correlation length L v and below by the viscous 
scale r). Passive scalar fluctuations generated at the scale 
L form increasingly finer structures due to velocity ad- 
vection and this process results in a net flux of scalar 
variance to small scales, where it is eventually smeared 
out by molecular diffusivity at a scale r<i- Here we will 
consider the case where these scales are ordered as fol- 
lows: L v ^> L ri,Td- In the range L r r<i the 
average scalar flux is constant and equals the average in- 
put rate: this is the well studied inertial-convective range 
where displays non-gaussian statistics and anomalous 
scaling [jj 0] . Conversely, at scales larger than L there 
is no scalar flux. Accordingly, one would expect Gaus- 
sian statistics and equipartition of scalar variance, i.e. 
the hallmarks of statistical equilibrium. This expecta- 
tion is correct at r > L„, where the dynamics of the 
large-scale passive scalar is ruled by an effective diffusion 
equation resulting in a Gaussian statistics. However, in 
the intermediate range L v > r > L, deviations from 
"thermal equilibrium" might arise as a consequence of 
turbulent transport. Indeed, Falkovich and Fouxon [j| 
have recently shown — in the context of the Kraichnan 
model of passive scalar advection where the velocity field 
is gaussian, self-similar and short-correlated in time Q 
— that the scalar field shows a highly nontrivial behav- 



ior at scales larger than the pumping correlation length 
L, with significant differences from Gibbs statistical en- 
semble. In this Letter we show that these results extend 
to passive scalar advection by a realistic flow, namely 
two-dimensional Navier-Stokes turbulence in the inverse 
cascade range. This flow has been studied in great de- 
tail both experimentally (in fast flowing soap films Q 
and in shallow layers of electromagnetical ly d riven elec- 
trolyte solutions and numerically [jj, [n| ■ The ve- 
locity is statistically homogeneous and isotropic, scale- 
invariant with exponent 1/3 (no intermittency correc- 
tions to Kolmogorov scaling) and with dynamical cor- 
relation times. This flow has also been utilized to inves- 
tigate passive scalar transport in the scalar flux range 
L > r > rd 0, and multi-particle dispersion, an 
intimately related subject From this ensemble 

of studies it emerged that the lessons drawn from the 
study of the Kraichnan model are indeed relevant for the 
qualitative understanding of passive scalar transport by 
realistic flows; and this will turn out to be true for the 
present case as well. 

Some basic, although incomplete, information about 
the scalar statistics in the supposed "thermal equilib- 
rium range" can be gained by studying the isotropic 
spectrum E(k) = 27rfc(|6>(fc)| 2 ) for kL > 1. Let us re- 
call that the equilibrium statistics for the scalar field 
would be described by the Gibbs functional V[9] = 
Z~ x exp[— (3 J \9(k)\ 2 dk], i.e. the Fourier modes should 
behave as independent Gaussian variables with equal 
variance 1/(2/3). As shown in Fig. ^ we observe E(k) ~ 
k, in agreement with equipartition arguments; moreover, 
the statistics of single Fourier modes is indistinguishable 
from Gaussian. However, we anticipate that from those 
findings alone one cannot state conclusively that large- 
scale passive scalar is in a thermal equilibrium state, 
given that they do not allow to rule out the possibility of 
long-range correlations for higher order observables (e.g. 
four-point scalar correlations). 

A more refined description of the large-scale properties 
of passive scalar can be obtained in terms of the coarse- 
grained field 



.(x,t) = J G r (x - y)0(y,t)dy 



(2) 



2 



where G r acts as a low-pass filter in Fourier space (for 
instance, the top-hat filter G r (x — y) = l/(irr 2 ) if 
\x — y\ < r and zero otherwise; or the Gaussian filter 
G r (x-y) = (2vrr 2 )exp[-|a;-y| 2 /(2r 2 )]). For r ->■ the 
filter reduces to a two-dimensional (5-function and there- 
fore^. — > 9. The statistics of 8 is typically supergaussian 
|l5l Il6| : its probability density function has exponential- 
like tails even for a gaussian driving force /. Indeed, in 
the latter case it can be shown that 9 is the product of 

two independent random variables 9 = FqT where 4> is 
a gaussian variable of zero mean and unit variance, Fq is 
the average injection rate of scalar fluctuations, and T is 
a positive-defined random variable, independent from 4>. 
The variable T is essentially the time taken by a spherical 
blob of minute initial size to disperse across a length L for 
a given flow configuration 0,0. For example, assum- 
ing a poissonian distribution for the rare events when T 
exceeds its average value yields exponential tails for the 
pdf of 9 0] . The distribution of 9 r is supergaussian as 
well; however, as r increases above the forcing correlation 
length, the probability density of 9 r tends to a gaussian 
distribution, as it is clearly seen by the scale-dependence 
of the distribution flatness and hyperflatness (see Fig.[5J). 

Within the framework of Gibbs statistical equilibrium, 
the scalar field has vanishingly small correlations above 
the scale L: therefore one could view 9 r as the sum of 
N Pi (r/L) 2 independent random variables (identically 
distributed as 9) divided by N. By central limit the- 
orem arguments |2(ij |. the moments of order 2n of the 
coarse-grained scalar field (odd-order moments vanish by 
symmetry) should then scale as N~ n , giving ((9 r ) 2n ) ~ 
(9 2 ) n (r I L)~ 2n . This is a very good estimate for n = 1: 
indeed, as shown in Fig. 03 the product (r/L) 2 (9 2 ) has a 
very neat plateau. This is consistent with the fast decay 
of the two-point scalar correlation (9(x, t)9(x + r, t)) at 
r > L. Indeed, in this case the second-order moment 
(9 2 ) = Jd yi dy 2 G r ( yi - x)G r (y 2 - x)(9( yi ,t)9(y 2 ,t)) 
is dominated by contributions with |yi — y 2 \ ^ L yield- 
ing (9 2 ) ~ (9 2 )(r / L)~ 2 . Alternatively, by Fourier trans- 
forming the coarse-grained field one obtains (Of) = 
f\G r (k)\ 2 (\9{k,t)\ 2 )dk ~ (9 2 )(r/L)~ 2 , since the trans- 
formed filter G r {k) is close to unity for kr -C 1 and falls 
off very rapidly for kr > 1, and \9(k : t)\ 2 ~ (9 2 )/{nL 2 ). 
In summary, two-point statistics appears to be consis- 
tent with Gibbs equilibrium ensemble. The situation for 
multi-point correlations will turn out to be different. 

A careful inspection of higher-order moments shows a 
less good agreement with central-limit theorem estimates 
(see Fig. |3Jl : this points to the existence of subleading 
contributions to the moments (9 2n ) for n > 1 arising from 
long-range correlations of multiple scalar products. In 
order to quantify more precisely the rate of convergence 
to gaussianity and its relationship to long-range correla- 
tions, it is useful to consider the cumulants of the random 
variable 9 r . According to central limit theorem [2(|, the 
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FIG. 1: (a) Passive scalar and velocity spectra. The data re- 
sult from the time integration of the two-dimensional Navier- 
Stokes equations dtv + v ■ Vw = — Vp + uAv + F — av and 
eq. Q by a pseudospectral method on a 1024 2 grid. The 
passive scalar is injected by a Gaussian, (^-correlated in time, 
statistically homogeneous and isotropic forcing restricted to a 
narrow band of wavenumbers. The initial condition for the ve- 
locity field is a configuration taken from a previous long-time 
integration and thus already at the statistically stationary 
state. The passive scalar starts from a zero field configura- 
tion, and after a transient of a few large-eddy turnover times 
t v — L v /v rma where L v is the integral scale of the velocity 
field, it reaches its own statistically steady state as well. Time 
averages are taken after this relaxation time has elapsed, for 
a total duration of more than 10 4 scalar correlation times 
tl ~ t v (L/ L v ) 2 ^ 3 . Here L/L v ~ 0.02. The velocity spectrum 
agrees with the Kolmogorov prediction k~ 5 ^ 3 and the passive 
scalar one follows very closely the equipartition spectrum in 
two- dimensions E(k) ~ k (see also the inset of panel (b)). (b) 
The marginal probability density function of a single Fourier 
amplitude 9(k) is indistinguishable from a Gaussian (dotted 
curve) for all wavenumbers in the range L _1 > k > L~ x . 
Here are shown three wavenumbers with kL — 0.5, 0.25, 0.12. 
In the inset is shown the spectral density (\6(k)\ 2 ) that shows 
a neat plateau at kL <C 1 (notice the linear scale on the ver- 
tical axis). 
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FIG. 2: The flatness and hyperflatness of the coarse-grained 
scalar field as a function of r/L, normalized by their gaussian 
values (2n — 1)!!. For r — > the curves tend to the flatness 
factors of the field 0: the numerical values correspond to a 
supergaussian probability density function \np{6) ~ — 8 1 ' 6 . 

cumulant of order 2n should vanish with jV~ 2n+1 leading 
to an expected scaling ((6* n )) ~ ((9 2n ))(r/ L)- 4n + 2 . Let 
us reiterate that the former expression is expected to be 
valid in absence of scalar correlations across lengthscales 
r>L. 

For n — 1 we have ((0 2 ,)) = {0 2 .) whose behavior 
has been already detailed above. In Fig. 01 we show 
the behavior of ((&*:)) = (9f.) - 3(9 2 ) 2 and ((0*)) = 
(6»f) - 15(6*2) {Of) + 30(6»2) 3 . For the fourth-order cumu- 
lant, we observe a scaling law very close to the theo- 
retical expectation ((fl*)) ~ ((0 4 ))(r/-L)- 16 / 3 @. This 
has to be contrasted with the scaling law (r/L)~ 6 given 
by central limit arguments. The breakdown of central 
limit theorem is due to the existence of long-range dy- 
namical correlations in the range r»L. These exclude 
the possibility of a true Gibbs statistical equilibrium at 
large scales. The leading contribution to the fourth-order 
cumulant ((&%)) = J dyi dy 2 dy 3 dy 4 G r (yi - x)G r {y 2 - 
x)G r {y 3 - x)G r {y A - x)((9( yi ,t)9(y 2 ,t)e(y 3 , t)9(y±, *)» 
comes from configurations with the four points arranged 
in two pairs of close particles (e.g. \yx — y 2 \ ^ L 
and 1 1/3 — 1/4 1 < L) separated by a distance r (e.g. 
1 2/i — j/3 1 ~ r). Otherwise stated, two-point correla- 
tors of the squared scalar field {{8 2 (x,t)9 2 (x + r,t))) 
must display a nontrivial scaling (r/L) -4 / 3 . We will 
get back to the issue of the statistics of 9 2 momentar- 
ily. The sixth-order cumulant ((9®)) is extremely difficult 
to measure because of the strong cancellations between 
various terms. Upon collecting the statistics over about 
ten thousand scalar correlation times, we can conclude 
that the results are consistent with the power-law decay 
((06)) ~ ((9 6 ))(r/L)- 22 / 3 suggested by the theory, and 
arising from terms like ((9f))((9r)) that appear in the 
expansion of the sixth-order cumulant 0. The actual 



FIG. 3: Moments of the coarse-grained scalar field (9^ n ) com- 
pensated by the thermal equilibrium expectation (r/L) 2n . 
The errorbars are determined by dividing the sample in ten 
subsamples and computing the dispersion around the mean. 
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FIG. 4: Cumulants of order 4 and 6 for the coarse-grained 
scalar field. The best fits for the slopes give exponents 5.32 ± 
0.05 for n = 2 and 7.5 ± 0.5 for n = 3. The theoretical values 
— 16/3 and —22/3 are shown for comparison. 

exponent for ((9f)} cannot be determined with great pre- 
cision, yet it lies within the range between —7 and —8, 
thus definitely different from the central-limit-theorem 
expectation —10. 

Further insight on the deviations from statistical equi- 
librium at large scales can be gained by studying the 
statistics of the coarse-grained squared scalar field 

0< 2 >(*,t) = / G r (x-y)9 2 (y,t)dy . (3) 

(2) 

The cumulants of 9 r give useful information about the 
presence of long-range correlations of the field 2 . The 

(2) (2} 

first-order cumulant ({$)■ )) = (#r ) is trivially equal to 
(0 2 ). The second-order cumulant ((6>f 2) )) = (oi 2) ) - 
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this is not the case for Ti and T2 because of the under- 
lying velocity field. Therefore, the long power-law tail 

for ((9 r 2] )) arises from events where (TiT 2 ) » (T) 2 - 
This amounts to say that two blobs of initial size smaller 
than L, released at a distance r > I in the same flow, 
do not spread considerably by turbulent diffusion (i.e. 
Ti, 2 > (T)) with a probability - (r/L)- 4 / 3 . 
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FIG. 5: Cumulants of the coarse-grained, squared scalar field 

(2) 

Of- . The definitions for the low-order cumulants of a generic 
random variable z are: ({z}} = (z), ({z 2 )} = (z 2 ) — (z) 2 , 
«z 3 » = (z s ) - 3(z)(z 2 ) + 2(z)\ ((z 4 )) = (z*) - i(z)(z 3 ) - 
3(z 2 ) 2 + 12(z 2 )(z) 2 -&(z) A . 




FIG. 6: A snapshot of the squared scalar field 6 2 . Remark 
the inhomogeneous distribution of scalar intensity originating 

from long-range correlations ((8^ )) ~ r _4//3 . 

(8 r 2 ^) 2 for a scalar field in thermal equilibrium should 
decay rapidly to zero at large scales r > L. On the con- 
trary, as shown in Fig. we observe a slow power-law 
decay with an exponent close to the theoretical expecta- 

tion«^ 2)2 ))c((0)(r/ J L)- 4 / 3 0. 

Higher-order cumulants behave self-similarly as 

((d (2)n )) ~ {{6 (2) )) n ~ 1 . This result can be interpreted in 
terms of the geometrical properties of the positive mea- 
sure defined by the squared scalar field: at scales r > L 
the field 9 2 appears as a purely fractal object with di- 
mension Dp w 2/3 (see Fig. |SJ). We end up by briefly 
discussing the physical origin of long-range scalar correla- 
tions. For a gaussian forcing we have 2 (xi, t)9 2 (x2, t) = 
Fq4>\4>\TiT 2 . At distances \x± — x 2 \ = r 3> L the two 
gaussian variables <j>\ and 4>2 are independent. However, 
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